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ABSTRACT 

We examine whether the ultraluminous infrared galaxies that contain a type I Seyfert nucleus (a type 
I ULIRG) are in the transition stage from ULIRGs to quasi-stellar objects (QSOs). To inspect this issue, 
we compare the black hole (BH) mass, the bulge luminosity and the far infrared luminosity among type I 
ULIRGs, QSOs and elliptical galaxies. As a result, we find the following results; (1) The type I ULIRGs 
have systematically smaller BH masses in spite of the comparable bulge luminosity relative to QSOs and 
elliptical galaxies. (2) The far-infrared luminosity of most type I ULIRGs is larger than the Eddington 
luminosity. We show that above results do not change significantly for 3 type I ULIRGs that we can 
estimate the visual extinction from the column density. Also, for all 8 type I ULIRGs, we investigate 
the effect of uncertainties of BH mass measurments and our sample bias, so that it turns out that our 
results do not alter even if we consider above two effects. In addition, Anabuki (2004) revealed that their 
X-ray properties are similar to those of the narrow line Seyfert 1 galaxies. These would indicate that 
active galactic nuclei (AGNs) with a high mass accretion rate exist in the type I ULIRGs. Based on all 
of these findings, we conclude that it would be a natural interpretation that type I ULIRGs are the early 
phase of BH growth, namely the missing link between ULIRGs and QSOs. Moreover, by comparing our 
results with a theoretical model of a coevolution scenario of a QSO BH and a galactic bulge, we show 
clearly that this explanation would be valid. 

Subject headings: galaxies: active — galaxies: bulges — galaxies:formation — galaxies:starburst — 
quasars: general — black hole 



1. INTRODUCTION 

Up to now, the Infrared Astronomical Satellite (IRAS) 
made the remarkable discovery of a new class of galax- 
ies, the ultraluminous infrared galaxies (ULIRGs, i.e., 
those having infrared luminosities greater than Lir(8 — 
1000/im) > I0 12 L Q ), which emit the bulk of their energy 
at infrared wavelengths. A lot of studies (e.g., Sanders & 
Mirabel 1996) have established that ULIRGs are gas-rich 
galaxies that in most cases have undergone a recent strong 
interaction with other galaxies eventually leading to a com- 
plete merger of the two. A commonly accepted explana- 
tion is that during such a galaxy collision the interstellar 
medium is transported toward the circumnuclear environ- 
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mcnt and is concentrated and compressed there, thus re- 
sulting in starburst on scales of less than a few kiloparsecs. 
In the next phase, the low angular momentum gas in the 
starburst region may fall into and accrete onto the central 
massive BH. In addition, the luminosity and space den- 
sity of ULIRGs are similar to those of QSOs. Moreover, 
the luminosity function of IRAS galaxies is of a double 
power-low type, and thus differs from that of the normal 
galaxies that show exponential rollover and is rather sim- 
ilar to that of QSOs or starbursts (Scoville 1992). Thus, 
it has been suggested that ULIRGs are powered by the 
heavily obscured QSOs (e.g., Sanders et al. 1988) or the 
starbursts (Joseph & Wright 1985). However, the physical 
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relation between ULIRGs and QSOs has been an issue of 
long standing. 

For the early type galaxies and QSOs, recent high 
spatial-resolution observations have suggested that the 
mass of a supermassive black hole (SMBH) tightly cor- 
relates with the mass, the velocity dispersion, and the 
luminosity of a galactic bulge (e.g., Kormendy & Rich- 
stone 1995; Richstone et al. 1998; Laor 1998; Tremaine et 
al. 2001; McLure & Dunlop 2001, 2002; Marconi & Hunt 
2003; Kawakatu & Umemura 2004). It has been found 
that the relatively low-z QSO hosts are mostly luminous 
and well-evolved early-type galaxies (e.g., McLeod & Rieke 
1995; Bahcall et al. 1997; Hooper, Impey, & Foltz 1997; 
McLoed, Rieke & Storrie-Lombardi 1999; Brotherton et al. 
1999; Kirhaokos et al. 1999; McLure et al. 1999; McLure, 
Dunlope, & Kukula 2000; Falomo et al. 2003; Dunlop 
et al. 2003). Recently, Veilleux, Sanders & Kim (1999) 
have shown that the percentage of AGNs is 30 — 50% for 
-Lir > 10 12 i Q . These findings suggest that the formation 
of a ULIRG, a SMBH, a galactic bulge and a QSO would 
be related to each other. However, this physical link is an 
open question. 

From the theoretical points of view, Kawakatu, 
Umemura & Mori (2003; hearafter KUM03) suggested a 
potential mechanism to build up a SMBH. They consider 
the effect of the radiation drag 6 , which extracts angular 
momentum from interstellar medium in starburst galaxies 
and thereby drives the mass accretion onto a galactic cen- 
ter (Umemura 2001; Kawakatu & Umemura 2002; Sato 
et al. 2004). On the basis of the radiation drag model, 
they proposed a new picture for a QSO formation. First, 
we regarded classical ULIRGs as the starburst galaxies, 
in which there is little AGN activity (if any). Next, we 
predicted the possibility of the "proto-QSO phase" , which 
is the optically thin and the total luminosity is dominated 
still by the bulge stars, although there is significant AGN 
activity in them. In this phase, a BH is still growing 
through the mass accretion and the BH-to-bulge mass ra- 
tio is smaller than that of QSOs and elliptical galaxies. 
And then, once a BH has grown fully and the central AGN 
dominates the total luminosity, the galaxy is regarded as 
a QSO (see also figures in KUM03). However, a proto- 
QSO has not been identified observationally yet although 
it is essential to clarify what objects correspond to proto- 
QSOs. 

Recently, ultra deep X-ray observations suggested that 
the submillimeter galaxies (SMGs) at z ~ 1 — 3 have the 
mass accretion rates approximately an order of magnitude 
lower than those of the coeval QSOs, assuming the Edding- 
ton luminosity (Alexander et al. 2005). Moreover, Borys 
et al. (2005) found that SMGs have smaller BH masses 
than QSOs with respect to the same mass range of bulges, 
and thus the SMGs may correspond to the "proto-QSO" 
predicted in KUM03. In the nearby galaxies, Canalizo & 
Stockton (2001) proposed that the infrared-selected type 
I AGNs (hereafter we refer them the "type I ULIRGs") 
are a transitional stage between ULIRGs and QSOs as 
their host galaxies are undergoing the tidal interactions or 
the mergers accompanied by the massive starbursts (see 
also Zheng et al. 2002; Lipari et al. 2005). In addition, 

6 The radiation drag in the solar system is known as the Poynting-Rc 
a similar effect on the formation of massive BHs (Umemura, Loeb, & 



most of them indicate the full width with half-maximum 
(FWHM) of the broad H/3 line less than 2000 km/s (Zheng 
et al. 2002), thus AGNs would be actually narrow line 
Seyfert 1 galaxies (NLSls). Recently, Anabuki (2004) 
studied 27 ULIRGs using the X-ray imaging and spec- 
troscopic observations with ASCA, Chandra, and XMM- 
Newton. Among their sample, 10 ULIRGs were identified 
as the type I ULIRGs. After correcting the absorption ef- 
fects, he found that seven luminous type I ULIRGs show 
the soft X-ray excess, the SED with a steep photon index 
(r2-iokcv > 2), and also the violent flux change (excess 
variance <~ 0.01 — 0.1), which are characteristic properties 
of NLSls. These X-ray properties would imply that AGNs 
with smaller BHs and high mass accretion rates exist in 
type I ULIRGs (e.g., Pound et al. 1995; Boiler et al. 1996; 
Mineshige et al. 2000). Moreover, Mathur et al. (2000) 
suggested that NLSls may be Syls in the early stage of 
their evolution if the BHs in NLSls are under massive with 
respect to their host bulges. Hence, it is likely that all of 
above characteristics support the hypothesis that the type 
I ULIRGs are the early phase of BH growth. However, the 
previous works have never examined that the BH-to-bulge 
relation among type I ULIRGs, QSOs and elliptical galax- 
ies. By investigating this issue, we can reveal whether type 
I ULIRGs have systematically smaller BH than QSOs and 
elliptical galaxies. By combining the BH-bulge relation for 
type I ULIRGs with the previous works for type I ULIRGs, 
we will test if they are really "proto-QSOs" , which are the 
transition phase from ULIRGs into QSOs. To this end, we 
demonstrate the relation among type I ULIRGs, QSOs, el- 
liptical galaxies on a BH mass (Mbh) versus an absolute 
i?-band magnitude of a galactic bulge (Mjt(bulge)). Also, 
we examine a BH mass (Mbh) and a far- infrared (FIR) 
luminosity (Lfir) for type I ULIRGs and QSOs, in order 
to constrain the origin of type I ULIRGs. 

The paper is organized as follows. In §2, we describe 
how type I ULIRGs, QSOs and elliptical galaxies are se- 
lected. In §3, we briefly review the technique of estimating 
BH masses from the broad emission-line widths for type 
I ULIRGs. In §4, we plot the data of type I ULIRGs in 
the Mbh — Mr (bulge) diagram and also compare them 
with that of QSOs and elliptical galaxies. Next, we show 
the Mbh — Lfik elation for type I ULIRGs and QSOs, 
and then compare with each others. Finally, we constrain 
the optical extinction of the central regions for type I 
ULIRGs by using the results of hard X-ray observations. 
In §5, we summarize our observational results for the type 
I ULIRGs, and then we compare them with one of the the- 
oretical models (KUM03 model). Section 6 is devoted to 
the conclusions. Through this paper, we adopt the Hub- 
ble parameter H =75 km s _1 Mpc -1 and the deceleration 
parameter go =0.5, and have converted the results from 
published papers to this cosmology to facilitate with com- 
parisons. 

2. SAMPLE SELECTION 

The aim of our study is to clear if type I ULIRGs are 
really the transition phase from ULIRGs into QSOs. To 
accomplish this, we need to use a type I ULIRG sample as 
the data of FWHM (H/3), the optical luminosity at 5100A 

son effect. Note that, in the early universe, Compton drag force has 
tier 1993). 
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in the rest frame and the R-band absolute magnitude of 
host bulges are coeval in these objects. For comparison, 
we also compile an optically-selected QSO sample and an 
elliptical galaxy sample, for which the BH mass and the 
i?-band absolute magnitude of a bulge are available. The 
details of these samples are given as follows. 

(1) The type I ULIRG sample is from Zheng et al. 
(2002). This sample was compiled from ULIRGs in the 
QDOT redshift survey (Lawrence et al. 1999), the Uy 
ULIRG survey (Kim & Sanders 1988), and an IR QSO 
sample selected from the cross-correlation of IRAS Point- 
Source Catalogue with the ROSAT All-Sky Survey cat- 
alog (Boiler et al. 1992). All the type I ULIRGs se- 
lected by Zheng et al. (2002) are ULIRGs with mid- 
infrared to far-infrared properties from IRAS observations. 
From this sample, we choose all 8 type I ULIRGs (IRAS 
F07599+6508, IRAS F11119+3257, IRAS Z11598-0112, 
IRAS F13342+3932, IRAS F15462-0450, IRAS F21219- 
1757, Mrk 231, and Mrk 1014), for which have both 
data of the width of broad H/3 line, the luminosity at 
5100A and the i?-band absolute magnitude of host galax- 
ies, Mr (bulge). Therefore, these objects are the best and 
maximal sample to achieve our aim at this time. As for 8 
type I ULIRGs, we obtain the FWHM (H/3) and the opti- 
cal luminosity at 5100A in the rest frame by Zheng et al. 
(2002) and M R (bulge) by Veilleux et al. (2002). In their 
paper, the contribution from the i?-band absolute magni- 
tude of AGNs were removed for Mr (bulge) (see Veilleux 
et al. 2002 for this procedure). In addition, according to 
Veilleux et al. 2002, all 8 type I ULIRGs are the single 
nucleus onjects. In Zheng sample, 23/25 type I ULIRGs 
have the optical luminosity at 5100A and the width of 
broad H/3 line. Then, we have compared the 8 selected 
objects with the rest 15 objects for these two properties. 
Figure 1 shows the optical luminosity at 5100A against the 
width of broad H/3 line. The filled red circles denote the 8 
selected type I ULIRGs, while the open red circles repre- 
sent the rest 15 type I ULIRGs whose bulge luminosities 
were not available. As seen in this figure, the 8 selected 
type I ULIRGs would be representative of large population 
with respect to the optical luminosity at 5100A and the 
width of broad H/3 line. In addition, the significant dif- 
ferences between our sample and the rest does not appear 
for the rages of the redshift (0.1 < z < 0.4) and infrared 
luminosity (Lir > 10 12 L Q ). 

(2) The optically-selected QSO sample comprises 29 
Palomar Green quasars (PG QSOs) from 30 luminous 
quasars (My < —23) published by McLure & Dun- 
lop (2001). This QSO sample consists of two optically 
matched subsamples of 17 radio-quiet QSOs and 13 radio- 
loud QSOs. The advantage of this sample is that all mem- 
bers have accurate bulge luminosities available from two- 
dimensional modeling of HST images. The average red- 
shift of the QSO sample is around 0.2. In this paper, 
we excluded PG 0157+001 (Mrk 1014) from the optically- 
selected QSO sample since it is categorized by type I 
ULIRGs. All 29 QSOs have the data of a BH mass and R- 
band magnitude of a bulge compiled by McLure & Dunlop 
(2001) and Dunlop et al. (2003). In their paper, the con- 
tribution from the i?-band absolute magnitude of AGNs 
were also removed for Mft(bulge). For 13 PG QSOs in our 
sample, their infrared flux was taken from Sanders et al. 



(1989) and Haas et al. (2000, 2003). 

(3) The elliptical galaxy sample consists of 20 objects 
drawn from the list of 37 nearby inactive galaxies with 
dynamical BH measurements published by Kormendy & 
Gebhardt (2002). In this paper, our main purpose is to in- 
vestigate the physical link between type I ULIRGs, QSOs 
and elliptical galaxies. Thus, we excluded those galaxies 
in the Kormendy & Dunlop list that were not E-type mor- 
phology (including lenticulars). The Kormendy & Geb- 
hardt list is made up of 20 E-type galaxies. All 20 ellipti- 
cal galaxies have the data of the BH mass and the 5-band 
absolute magnitude of the bulge (Kormendy & Gebhardt 
2002; Gebhardt 2003). To convert the B-band magnitude 
to the J?-band, standard bulge colors of B-R=1.57 were 
assumed (Fukugita et al. 1995). 

We summarized the various physical parameters of type 
I ULIRGs, QSOs and elliptical galaxies in Table 1. 
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Fig. 1. — Optical luminosity at 5100A against the width of broad 
H/3 line. The filled red circles denote the 8 selected type I ULIRGs, 
while the open red circles represent the rest 15 type I ULIRGs. The 
distribution of the 8 selected type I ULIRGs are similar to that of 
the rest. 

3. BLACK HOLE ESTIMATE OF TYPE I ULIRGS 

As mentioned §2, the BH masses of QSOs and elliptical 
galaxies have been already estimated by previous works, 
but we do not have those of 8 type I ULIRGs. Thus, 
we need to estimate them within the present paper. The 
method to estimate a BH mass is based on the assump- 
tion that the motion of ionized gas clouds moving around 
the black hole is dominated by the gravitational force and 
the clouds within the broad line region (BLR) is virialized 
(e.g., Peterson & Wandel 1999, 2000). Thus, the central 
black hole mass can be expressed by Mbh ~ Rblrv 2 /G, 
where v is the velocity dispersion of matter at the size of 
the broad line region -Rblr, which is gravitationally bound 
to the BH. Then, the central mass can be estimated as 



M BH = 1-5 x 10 b 



-Rblr 



f- 



"FWHM 



M . (1) 



It — days j \ 10 3 km s 
The velocity dispersion v can be estimated from the 
FWHM of H/3 broad line emission v = /wfwhm, by as- 
suming the BLR gas is in isotropic motions (/ = \/3/2). 
Based on 17 Seyfert galaxies and 17 optically-selected PG 
QSOs, Kaspi et al. (2000) determine an empirical rela- 
tionship between the size of the broad- line region, -Rblr, 
and optical continuum luminosity, ALA(5100A) rcst , where 
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-Rblr is the distance of the emission-line clouds respond- 
ing to the central continuum variation as determined from 
reverberation mappings (see Kaspi et al. 2000): 



R 



BLR 



(32.9 



+2.0\ 
1.9,7 



AL A (5100A) 

r 



10 44 ergs" 



0.70±0.033 



light— days. 
(2) 

By combining the equation (1) and (2), we can obtain 
the following formula: 



M BH = (4.9t£S) x 10 



AL A (5100A) rost 



0.70±0.033 



10 44 ergs 

/ fFWHM 



T ) 2 M Q . (3) 



V 10 3 kms 

As seen from the equation (3), the luminosity at 5100A in 
the rest frame and the FWHM of H/3 are needed to esti- 
mate a BH mass. We should also keep in mind that the 
equation (2) holds not only for broad-line type I AGNs, 
but also for NLSls (Peterson et al. 2000). In order 
to evaluate the optical luminosity at 5100A in the rest 
frame, L A (5100A) rost , we use the formula LA(5100A) rost = 
471x^(1 + z)F\ (5100(1 + z)A) obs , where c?l is the luminos- 
ity distance which is given by c?l = (cz/H ){l + z/A) for a 
small value of z. Here, the FWHM (H/3) and the observed 
flux at 5100A, F A (5100(1 + z)A) obs , are given by Zheng et 
al. (2002) and are measured directly from their spectra. 

In this paper, we assume that the contribution from 
central AGNs dominate the optical emission at 5100A 
in the rest frame and that the stellar continuum emis- 
sion is negligible. To confirm if this assumption is rea- 
sonable, we check the flux ratio of the optical emis- 
sion at 5100(1 + z)k and the hard X-ray emission of 2- 
10(l+z)keV for 5 type I ULIRGs. We get the flux of 
hard X-ray from Anabuki (2004) . As a consequence, 3 ob- 
jects (IRAS F11119+3257, IRAS Z11598-0112, and Mrk 
1014) have nearly same flux ratio as that of PG QSOs, 

10 3 ~ 5 . As for PG QSOs, we 

the data of i<2-io(i+z)keV 
2000). 

Although other 2 objects (IRAS F07599+6508 and Mrk 
231), which are the broad absorption line QSOs (BAL 
QSOs), have an extremely low luminosity at the X-ray 
band, it has been considered that it is not be due to the 
intrinsic effect but the absorption effect (e.g., Gallagher et 
al. 2002). Thus, it is expected that these 2 BAL QSOs 
have almost same flux ratio as that of PG QSOs intrin- 
sically. Since the contributions from AGNs dominate the 
optical emission for PG QSOs, it is valid assumption that 
the central AGNs in type I ULIRGs dominantly power the 
optical emission at 5100A in the rest frame. All the basic 
parameters of 8 type I ULIRGs are listed in Table 2. 

4. RESULTS 

4.1. M B h - M R Relation 

Figure 2 plots the i?-band absolute magnitude of bulge 
components (spheroidal components), Mr (bulge) [mag] , 
versus the black hole mass, Mbh [Mq] for 8 type I 
ULIRGs, 29 QSOs and 20 elliptical galaxies. In this fig- 
ure, the red circles represent type I ULIRGs, the blue 
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obtain this flux ratio from 

(George et al. 2000) and F 5100 ^ 1+ ^ (Kaspi et al. 



squares show QSOs, and the green circles denote ellip- 
tical galaxies. The solid line is the best-fitting relation 
for QSOs. This relation is given by log(MBH/M ) = 
-0.61(±0.08)M R (mag) - 5.47(±1.82), which is compara- 
ble to that found by Laor (1998). The underlines below the 
name of objects denote that they have the similar prop- 
erties at the soft and the hard X-ray band as those of 
NLSls (Anabuki 2004). As for all QSOs and ellipticals, 
they corrected the effect of Galactic and internal extinc- 
tion (McLure & Dunlop 2001; Dunlop et al. 2003; Geb- 
hardt et al. 2003; Kormendy & Gebhardt 2003). In the 
case of type I ULIRGs, the effect of the visual extinction 
may not be negligible because of significantly larger FIR 
luminosity, which would imply the plenty of dusty gas. 
Thus, the optical extinction may affect on the estimation 
of the BH masses and the bulge luminosities. Then, we 
exhibit the extinction effect toward the BLRs of type I 
ULIRGs as the arrows in Figure 2. Hereafter, we use Ay 
as the total extinction toward the BLRs of type I ULIRGs. 
Also, the effect of the optical extinction for the hosts of 
type I ULIRGs may be also significant (e.g., Veilleux et 
al. 2002). If this is the case, then we would underestimate 
the bulge luminosity of type I ULIRGs. As for the mor- 
phology of host galaxies, a lot of works have showed that 
the surface brightness profiles of QSO hosts resemble the 
de Vaucouleurs profile (e.g., McLure et al. 1999). As for 
8 type I ULIRGs, the surface brightness profiles of 6 host 
galaxies in type I ULIRGs can fit with the de Vaucouleurs 
profile (elliptical-like) for R and K'-b&nd, and the others 
(IRAS 11598 and IRAS F15462) can fit with both the de 
Vaucouleurs profile and exponential profile (disk-like) for 
their band (Veilleux, Kim & Sanders 2002). Thus, we put 
on the label "(E/D)" for the latter 2 objects in Figure 
2. In this paper, we focus on their absolute magnitude 
of spheroidal components. If the magnitude of galaxies 
are dominated by disk components, our estimations of the 
i?-band absolute magnitude in bulges would be overesti- 
mated. Note that we assume / = V3/2 to estimate the 
BH mass for QSOs, while Mclure & Dunlop used / = 3/2, 
and thus our estimation for QSO BHs is three times as 
small as their estimation. 

As seen in Figure 2, we have found that type I ULIRGs 
have systematically smaller BH mass than QSOs and el- 
liptical galaxies in spite of the comparable bulge luminos- 
ity, if the visual extinction effect is small (Ay < 3) for 
type I ULIRGs. Namely, the BH mass ranges are Mbh ~ 
10 6 - 8 M Q for type I ULIRGs and M BB w 10 8 ~ 9 Af Q for 
QSOs and elliptical galaxies. However, if Ay of all type I 
ULIRGs is larger than ~ 3, the BH mass of type I ULIRGs 
would be similar to that of QSOs. Thus, in order to justify 
if a BH mass of type I ULIRGs is systematically small, we 
need to constrain the effect of the visual extinction, which 
will be discussed in §4.3. On the other hands, if the extinc- 
tions for their host galaxies are significant, then the data 
points just move on top in Figure 2. In short, this effect 
makes the difference between the BH-mass distributions 
of type I ULIRGs, QSOs and ellipticals large. In addition, 
we find that the elliptical galaxies are located at slightly 
lower parts than QSOs in Figure 2. According to KUM03, 
this result may indicate that the hosts of QSOs are slightly 
younger than that of elliptical galaxies (for details, see §5). 
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Fig. 2. — Absolute i?-band bulge magnitude versus black hole 
mass for 8 type I ULIRGs (red circles), 29 QSOs (blue squares) and 
20 elliptical galaxies (green circles). The black hole masses for type 
I ULIRGs are derived from their broad H/3 line widths by using 
equation (3). The black hole masses for QSOs are given by McLure 
& Dunlop (2001) and Dunlop et al. (2003). The black hole masses 
for elliptical galaxies are the dynamical estimates as Kormendy & 
Gebhardt (2003) and Gebhardt et al (2003). The black horizontal 
arrow shows the extinction effects for the BLR of type I ULIRGs. 
The under line below the name of objects denote that they have the 
same properties at hard X-ray band like NLSls (Anabuki 2004). The 
solid line is the best-fitting relation for the optical QSOs, which is 
log(M BH /M Q ) = -0.61(±0.08)M R (mag) - 5.47(±1.82). The sym- 
bols "(E/D)" denote the objects can fit with both de Vaucouleurs 
profile (elliptical-like) and exponential profile (disk-like) for R and 
A" band (Veilleux et al. 2002). 



4.2. A/bh-£fir Relation 

Figure 3 shows that the BH mass-to-FIR luminosity re- 
lation, A^bh — ^fir[40— 500^m] using the derived BH mass 
(see Table 1) and the FIR luminosity for 8 type I ULIRGs 
(Zheng et al. 2002). In order to compare the result of type 
I ULIRGs with that of QSOs, we select 13 PG QSOs from 
the sample selected by McLure & Dunlop (2001) and Dun- 
lop et al. (2003). Ah 13 PG QSOs have the data of the 
IRAS flux densities at 60 and 100/im (Sanders et al. 1989; 
Haas et al. 2000, 2003). As for 13 PG QSOs, we calculated 
their far-infrared luminosities following formula (Sanders 
& Mirabel 1996) , based on the flux densities from the IRAS 
Faint Source Catalog: Lfir[40 — 500yum] = 47rdi J Ci ; FTR, 
where the scale factor C(= 1.4— 1.8) is the correction fac- 
tor required to account principally for extrapolated flux 
longward of the IRAS 100 fim filter, and Ffir is defined 
as 1.26 x 10- 14 x (2.58/ 6 o + /ioo)[Wm^ 2 ] with / 60 and / 100 
being the IRAS flux densities at 60 and 100 /jm in unit of 
Jy. Here, we employ C = 1.8. The red circles show the 
type I ULIRGs and the blue squares and arrows represent 
QSOs. The black horizontal arrow shows the optical ex- 
tinction effect for the BLRs of type I ULIRGs. The thick 
solid line denotes the luminosity ratio of far infrared to 
Eddington luminosity (ipm/LEdd) equals unity. The thin 
solid, dashed and dot-dashed line are iFiR/^Edd = 10, 
Lpm/LEdd = 0.1 and L Fm /L Edd = 0.01, respectively. 



Finally, we check the effects of systematic errors in the 
methods used to evaluate BH masses and the effect of our 
sample bias. First, we discuss the systematic errors in 
methods. As for the type I ULIRGs, the uncertainties of 
the BH mass were estimated by error propagation using 
the optical luminosity at 5100A and FWHM of H/3 mea- 
surements given by Zheng et al. (2002) . The mean error of 
the BH mass is a factor 1.3. Although McLure & Dunlop 
(2001) did not show the uncertainties of BH masses clearly, 
in general the BH mass in the this way (see §3) is accurate 
within a factor 2-3 (e.g., Wang & Lu 2001; Marziani et al. 
2003; Shemmer et al. 2004). As seen in Figure 2, except 
for Mrk231, the BH mass of type I ULIRGs are ten times 
as small as that of QSOs and elliptical galaxies at the fixed 
i?-band magnitude of bulge. Furthermore, the mean error 
of BH mass for elliptical galaxies is a factor 3. Thus, it 
would be clear that the systematic difference of BH masses 
between type I ULIRGs and QSOs (or ellipticals) can not 
be explained by only systematic errors in methods of BH 
mass measurements. On the other hands, as shown in §2, 
the 8 selected type I ULIRGs would be representative of 
large population. Thus, the rage of BH masses in other 12 
type I ULIRGs are similar to that of the 8 selected type 
I ULIRGs with A/ B h = 10 6 ~ 8 Af Q . The rest three objects 
have massive BHs with ss 10 9 M Q . As a consequence, 12 
type I ULIRGs would be located farthest from the loca- 
tion of QSOs if the range of i?-band magnitude of bulge 
is from -22 to -24 that is typical range for the 8 selected 
type I ULIRGs. Thus, we found that the systematic dif- 
ference we found in Mr (bulge) — Mbh diagram would not 
be the effect of our sample bias. We should keep in mind 
that othee results we will show later (Figure 2 and 3) also 
do not change significantly by the systematic errors in the 
BH measurements and the effect of our sample bias. 



£ Mrk231 
07599 




0.1L Edd / 

./'O.OlZea 



Log(M BH /M n ) 



Fig. 3. — Far-infrared luminosity versus the black hole mass 
for 8 type I ULIRGs (red circles) and 13 PGQSOs (blue squares 
and arrows). The black horizontal arrow shows the effect of visual 
extinction for the BLRs of type I ULIRGs. The thick solid line de- 
notes the luminosity ratio of far-infrared to Eddington luminosity 
(iriR/^Edd) equals unity. The thin solid, dashed and dot-dashed 
line are LFiR/L E dd = 10, ipiR,/iEdd = 0.1 and LpiR/L E dd = 0.01, 
respectively. 

As seen in Figure 3, it turns out that the FIR luminos- 
ity is larger than the Eddington luminosity for most type 
I ULIRGs in Figure 3. By contrast, the FIR luminosity is 
more than one order of magnitude smaller than Eddington 
luminosity for QSOs, namely Lfir < O.lLEdd- This may 
indicate AGNs in type I ULIRGs with a high mass ac- 
cretion rate, or the existence of another power source, for 
which the promising candidate is a star burst. On this is- 
sue, we will discuss §5.1. In addition, Figure 3 shows that 
type I ULIRGs have systematically a smaller BH mass 
than QSOs at fixed the FIR luminosity. On the other 
hands, the large far- infrared luminosity would imply the 
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plenty of dusty gas (e.g., Haas et al. 2003). Hence, for type 
I ULIRGs it would show that the mass ratios of the total 
dusty gas to the BH arc much larger than that of QSOs. 
It would indicate that BH masses of type I ULIRGs have 
growth potential. 

4.3. Optical Extinction Effect 

As mentioned in §4.1 and §4.2, we found that type I 
ULIRGs have systematically smaller BHs than that of 
QSOs at given the absolute magnitude of bulge and given 
the FIR luminosity. However, the BH masses of type I 
ULIRGs would depend on the optical extinction effect, 
namely, the heavily extinction makes the derived BH mass 
smaller (see the horizontal arrow in Figure 2 and 3). On 
the other hands, it is well-known that Balmer decrement 
is a standard evaluation of optical extinction of the nar- 
row line regions (e.g., Osterbrock 1989). However, this 
method is generally invalid to estimate the amount of the 
visual extinction toward BLRs, because the flux ratios of 
broad Balmer lines are sometimes seriously affected by col- 
lisional excitation effects. As a complementary approach, 
the hard X-ray observations enable us to make a measure- 
ment of the extinction toward the nucleus, and to con- 
strain the BH mass more reliably. We should note that 
many observations suggested that the absorption column 
density (iVu) derived from hard X-ray is systematically 
large relative to optical extinction of Ay under the as- 
sumption of Galactic gas/dust mass ratio (e.g., Maiolino 
et al. 2001a; Watanabe et al. 2004). Thus, we esti- 
mate the optical extinction of 3 type I ULIRGs (IRAS 
F11119+3257, IRAS Z11598-0112 and Mrkl014) with the 
relation Ay/N n = 4.8^ 4 6 1 x 10" 23 magc m - 2 (la dis- 
persion), which is derived from the AGN observations 
(Maiolino et al. 2001b). As for two BAL QSOs (IRAS 
F07599+6508 and Mrk 231), the Ay/N B relation for nor- 
mal QSOs would not hold on BAL QSOs. However, re- 
cent works have found that the optical extinction of BAL 
QSOs is around 0.1-1 by comparing the composite non- 
BAL QSO spectra with the composite BAL QSO spectra 
(e.g., Brotherton et al. 2001; Richard et al. 2003). There- 
fore, the optical extinction of two BAL QSO in our sample 
may be less than Ay — 1. We summarize the visual ex- 
tinctions for type I ULIRGs in Table 3. As for at least 3 
type I ULIRGs to be able to estimate the visual extinction 
from the column density, we show that our results (§4.1 
and §4.2) do not change drastically. 

4.4. Summary of observational results 

Based on the above results, we summarize our findings 
on type I ULIRGs as follows. (1) The BH mass of type I 
ULIRGs is systematically smaller than QSOs and ellipti- 
cal galaxies despite of the comparable bulge luminosity to 
them. (2) Most type I ULIRGs have particularly large FIR 
luminosity against the Eddington luminosity. We show 
that above results do not change significantly for 3 type 
I ULIRGs even if we consider the effects of the visual ex- 
tinction. Also, for 8 type I ULIRGs, we investigate the 
effect of uncertainty of BH mass measurements and our 
sample bias, so that we found that our results do not al- 
ter. Additionally, their X-ray luminosity properties are 
similar to those of NLSls, whose X-ray properties reflect a 
high mass accretion rate (Anabuki 2004 in details). From 



all these findings, it would be a natural explanation that 
type I ULIRGs are the early phase of BH growth, namely 
the transition stage from ULIRGs to QSOs. In the next 
section, we will investigate whether this interpretation is 
reasonable by comparing our results with the theoretical 
prediction. 

5. DISCUSSIONS 

Based on our observational results for 8 type I ULIRGs, 
we have suggested that type I ULIRGs would be the tran- 
sition stage from ULIRGs to QSOs. However, it has not 
been cleared that the physical relationship between type 
I ULIRGs, QSOs and elliptical galaxies with the objec- 
tive of the coevolution of host galaxies and SMBHs. For 
this purpose, we compare our results with the theoretical 
predictions of a coevolution scenario of the galactic bulges 
and SMBHs. 

5.1. A Coevolution Scheme for SMBHs and Galactic 

Bulges 

Recently, KUM03 have constructed a physical model for 
a coevolution of a QSO BH and an early-type host galaxy. 
This model is based on the radiation drag incorporating 
the realistic chemical evolution that reproduces the color- 
magnitude relation of present-day bulge. Here, we briefly 
review the essence of their model as a pre-arrangement 
of the following discussions. In their model, they used 
an evolutionary spectral synthesis code 'PEGASE' (Fioc 
& Rocca-Volmerange 1997), in order to treat the realis- 
tic chemical evolution of the galactic bulges. According 
to results of KUM03, after a galactic wind epoch t w , the 
bolometric luminosity is shifted from the host-dominant 
phase to the AGN-dominant phase (the QSO phase) at 
the transition time i C rit- The former phase (t w < t < i C rit) 
corresponds to the early stage of a growing BH, because 
the mass accretion rate during this phase is so high that 
the mass growth of BH is significant. They defined this 
phase as a "proto-QSO". The proto-QSO phase is pro- 
ceeded by an optically thick phase before the galactic 
wind, which would correspond to a classical ULIRG. In 
this phase, they predicted that the BH is much smaller 
than the QSO phase. After the AGN luminosity exhibits 
a peak at i C ross, it fades out abruptly because almost all of 
the matter around BH has fallen onto the central BH. The 
fading nucleus could be a low luminosity AGN (LLAGN). 

5.2. Are Type I ULIRGs Missing Link Between ULIRGs 

and QSOs ? 

By using KUM03 model, we predict the evolution of 
Mr (bulge) and -Mbh for the different masses of bulges 
in Figure 4. Here, we calculate the evolutional tracks 
in the Mr-Mbh diagram for 4 different masses of bulges 
(M g0 = 10 10 , 10 11 , 10 12 , and2 x 1O 12 M ), where M g0 is 
the initial gas mass in galactic bulges. In this figure, we 
assume a Salpeter-type initial mass function (IMF) to be 
<j> = dn/dlogm* = A(m*/M ) -1 ' 35 for a mass range of 
[O.IMq, 6OM ]. The star formation rate (SFR) per unit 
mass at time t, C(t), is assumed to be proportional to the 
gas mass fraction / g , C(f) = kf g at t < t w and at t > t w 
C(t) = 0, where k is the constant rate coefficient. In Table 
4, we summarize the model parameters (M g0 , k, t w , i cr i t , 
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and tcross)- 7 The evolutions of different stellar masses 
of bulges proceed from left (smaller BH) to right (larger 
BH) in Figure 4. The black horizontal arrow denotes the 
optical extinction effect for the BLR of type I ULIRGs. 




Log(M BH /M Q ) 

Fig. 4. — Absolute ij-band bulge magnitude versus the black hole 
mass for 8 type I ULIRGs (red circles), 27 QSOs (blue squares), 2 
"unusual" QSOs (magenta squares) and 20 elliptical galaxies (green 
circles). The dashed arrow denotes the BH mass if we use the broad 
line width with 8000 km/s. The black horizontal arrow shows the op- 
tical extinction for the BLR of type I ULIRGs. The grids represent 
the prediction based on the radiation drag model for different mass 
of bulges (M g0 = 10 10 , 10 11 , 10 12 , and2 X 10 12 M Q ), where M g0 is 
the initial gas mass in galactic bulges. Each evolutions shift from 
left (smaller BH) to right (larger BH) in this figure. i w is the galac- 
tic wind time scale. f cr it is the time when the luminosity of bulges is 
equal to that of AGNs. tcross is defined as the time when almost all 
of the matter around BH has fallen onto the central BH. to denotes 
the final stage of the galaxy evolution. The host luminosity domi- 
nant phase (proto-QSOs) correspond to the region (t w < t < t cr it)- 
The AGN luminosity dominant phase (QSOs) correspond to the 
area (i C rit < t < tcross). We call the phase (tcross < t < to) a l° w 
luminosity AGN (LLAGN). 

As seen in Figure 4, it is found that most type I ULIRGs 
are located near the proto-QSO phase. It would indicate 
that type I ULIRGs are the early phase of BH growth 
within younger bulge and then they evolve into QSOs. 
Namely, it would suggest that type I ULIRGs are the 
missing link between ULIRGs and QSOs. In addition, we 
also find that two "unusual" QSOs (PG 0204+292 and 
PG 2247+140) are located near the region of the proto- 
QSOs. In fact, these "unusual" QSOs have a narrower H/3 
line width (< 2500km/s) among the samples of PG QSOs 
(McLure & Dunlop 2001). Thus, these may correspond 
to the objects on the way of evolving into QSOs. How- 
ever, as for the PG 0204+140, we should comment that 
it has not only a strong narrow component, but also a 
very broad line width (about 8,000 km/s: see Figure 6 in 
McLure & Dunlop 2001). If we estimate its BH mass by 
using larger FWHM, PG 0204+140 has the BH mass with 
3 x 10 8 Mq (the dashed arrow). If this is case, it is located 
near the region of normal QSOs in Figure 4. Additionally, 
our model can explain the difference from QSOs and ellip- 
tical galaxies as the evolution of host galaxies with nearly 
same BH mass. Moreover, according to KUM03, it would 
be expected that parent galaxies of type I ULIRGs are 
the starburst galaxies, in which there is any AGN activity. 
This corresponds classical ULIRGs. In this phase, their 



BH mass would be much smaller than the BH mass of type 
I ULIRGs. To sum up, we found that classical ULIRGs — 
type I ULIRGs — QSOs — ellipticals would be explained by 
the evolution sequence of the spheroidal systems. How- 
ever, these are some discrepancy between the predictions 
of KUM03 model and our results for the type I ULIRGs. In 
KUM03, they predicted that a proto-QSO is an optically- 
thin and the total luminosity is still dominated by the host 
components. On the other hands, as for type I ULIRGs 
we revealed the followings; (1) The AGNs powers the op- 
tical emission rather than the stellar components (see §3). 
It might imply that AGNs contribute the energy source 
of FIR luminosity. (2) Larger FIR luminosity than the 
Eddington luminosity indicates the optical depth of their 
hosts is still thick (§4.2). It is useful to consider the rea- 
sons of two disagreements. Hence, we will discuss in next 
sections. 

5.2.1. Origin of the FIR luminosity 

In KUM03 model, they assume the AGN luminosity to 
be the Eddington one. But, recent optical and X-ray ob- 
servations have suggested that almost of all type I ULIRGs 
would be the high mass accretion phases (Anabuki 2004; 
Hao et al. 2005). If the accretion rate onto a BH is 
the super-Edddington (M > Me), then the solution of 
accretion is optically-thick and called a slim-disk (e.g., 
Abramowicz et al. 1998), where the Eddinton accretion 
rate M E = L E /c 2 (e.g., Kato, Mineshige & Fukue 1998). 
In this type of disk, the photon-trapping effect (e.g., Katz 
1977; Begelman 1978) in the accretion flow plays an impor- 
tant role. In fact, the BH accretion luminosity can achieve 
up to w 10L E (Ohsuga et al. 2002, 2003, and 2005). 
Hence, by considering such a super-Eddington accretion 
in the type I ULIRGs, it is expected that the AGN lumi- 
nosity become larger than the host luminosity, and then 
this may be consistent with the observational result. If this 
is the case, moreover, then it turns out that the FIR lumi- 
nosity of type I ULIRGs can be explained by only the BH 
accretion luminosity (see thin solid line ipir/iEdd = 10 in 
Figure 3). On the other hands, the recent high resolution 
multiwavelength observations have indicated that the vast 
majority of ULIRGs are strongly interacting or merging 
galaxies (e.g., Clements et al. 1996; Murphy et al. 1996; 
Veilleux, Kim, & Sanders 2002). Thus, such a interaction 
or merger may cause a starburst (e.g., Larson & Tinsley 
1978; Noguchi 1988; Mihos & Hernquist 1994; Sanders 
& Mirabel 1996), and thus a starburst phenomenon may 
contribute to the FIR luminosity of the type I ULIRGs. 
However, there's plenty of scope for discussion whether the 
interactions or mergers trigger the starbursts as Bergvall 
et al. (2003) claimed that the interactions and mergers do 
not trigger the starburst. Nevertheless, it is worth empha- 
sizing again that the AGN luminosity only can produce 
the FIR luminosity of type I ULIRGs if the accretion rate 
is M > Me- This would impact on one of the important 
issues for ULIRGs, that is, whether the dominant energy 
sources of ULIRGs are the dust-obscured AGNs or star- 
bursts. 

5.2.2. Obscuring problem 



7 As for the effect of star formation history, if the mass range and slope of IMF and the SFR are changed to satisfy the spectrophotomctric 
properties of galactic bulges, the final BH mass is altered by a factor of ±50% (Kawakatu & Umemura 2004). 
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KUM03 was based on the monolithic model for a co- 
evolution of a SMBH and a galaxy. As we mentioned in 
§5.2.1., however, almost all ULIRGs have undergone the 
mergers or the interactions, and also a AGN phenomenon 
appears at the final merging stage (e.g., Clements et al. 
1996; Kim, Vcillcux, & Sanders 1998; Wu et al. 1998; 
Zheng et al. 1999; Canalizo & Stockton 2001; Cui et al. 
2001) Thus, the interactions or mergers would play a sig- 
nificant role in order to elucidate the physical link between 
ULIRGs and QSOs. If dusty galaxies interact or merger 
with optically thin proto-QSOs, then the dusty gas may 
be winded up and then the dust-enshrouded AGNs may 
form. The objects may observe as type I ULIRGs. After 
the dusty gas would be swept away by the multiple super- 
novae and AGN feedbacks (e.g., Granato et al. 2004), type 
I ULIRGs would evolve into the optically thin QSOs. As 
the other possibility, we may consider that type I ULIRGs 
may be the on the way to blowing away the dusty gas 
even within the monolithic scenario. Indeed, the duration 
of sweeping away the dusty gas is finite, so that the optical 
depth of bulges would be mildly optically thick during this 
phase, though we supposed that the galactic wind sweeps 
the dusty gas in an instant in Figure 4. Thus, in order to 
understand the physical properties of type I ULIRGs in de- 
tails, we will elucidate the wind effects and interactions or 
mergers with the dusty galaxies by using the sophisticated 
simulations in the future works. 

6. CONCLUSIONS 

We have investigated whether type I ULIRGs are really 
in the transition stage from ULIRGs to QSOs, using the 
data of infrared, optical and X-ray observations. To reveal 
this issue, we compare the BH mass, the bulge luminosity 
and the FIR luminosity among type I ULIRGs, QSOs and 
elliptical galaxies. Our main conclusions are: 

1 . The type I ULIRGs have systematically smaller BH 
masses in spite of the comparable bulge luminosity 
relative to QSOs. 

2. The far-infrared luminosity of most type I ULIRGs 
is larger than the Eddington luminosity. 



3. We have shown that our results do not change sig- 
nificantly for 3 type I ULIRGs that we can evalu- 
ate the effects of the visual extinction. In addition, 
we have examined the effect of uncertainty of BH 
mass measurements and our sample bias, so that 
we found that our results do not alter drastically. 

4. From the X-ray properties of type I ULIRGs, which 
are similar to those of the narrow line Seyfert 1 
galaxies (Anabuki 2004), these would indicate that 
AGNs with a high mass accretion rate exist in the 
type I ULIRGs. 

5. Based on all of these findings (1-4), it would be a 
natural interpretation that the type I ULIRGs are 
the early phase of BH growth, namely the miss- 
ing link between ULIRGs and QSOs. Moreover, by 
comparing our results with a theoretical model of 
a coevolution scenario of a QSO BH and a galactic 
bulge, we have found that this explanation would 
be reasonable. 

6. The AGN luminosity only can produce the FIR lu- 
minosity of type I ULIRGs if we consider the super- 
Eddington accretion. This result would impact on 
the origin of the infrared luminosity in ULIRGs. 
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Table 1. The Variable Physical Parameters 



Name 


log(JW B H/M ) 


Mr (host) 


log(L F m/-L©) 


references 


(i) 


(2) 


(3) 


(4) 


(5) 


Typo I ULIRGs 


IRAS F01572+0009 (Mrk 1014) 


7.78 


-23.76 


12.30 


1.2 


IRAS F07599+6508 


8.08 


-24.68 


12.10 


1,2 


IRAS F11119+3257 


7.11 


-22.61 


12.29 


1,2 


IRAS Z11598-0112 


6.28 


-22.61 


11.56 


1,2 


IRAS F12540+5708 (Mrk 231) 


7.92 


-22.49 


12.25 


1,2 


IRAS F13342+3932 


6.49 


-22.87 


12.14 


1,2 


IRAS F15462-0450 


6.26 


-20.91 


12.00 


1,2 


IRAS F21219-1757 


7.08 


-22.11 


11.67 


1,2 


PG QSOs 


PG 0137+012 


8.58 


-23.16 




3 


PG 0736+017 


8.01 


-22.70 




3 


PG 1004+130 


9.11 


-23.22 


< 11.76 


3,7 


PG 1020-103 


8.37 


-22.48 




3 


PG 1217+023 


8.42 


-22.83 




3 


PG 1226+023 


8.62 


-23.52 


12.12 


4,5 


PG 1302-102 


8.31 


-22.62 


12.09 


4,6 


PG 1545+210 


8.94 


-22.42 


< 11.38 


4,7 


PG 2135-147 


8.95 


-22.62 




3 


PG 2141 + 175 


8.75 


-22.92 




3 


PG 2247+140 


7.60 


-22.92 




3 


PG 2349-014 


8.79 


-23.52 


< 11.57 


3,7 


PG 2355-082 


8.40 


-22.72 




3 


PG 0052+251 


8.29 


-22.12 




3 


PG 0054+144 


8.91 


-22.72 




3 


PG 0204+292 


6.68 


-22.42 




3 


PG 0205+024 


7.86 


-20.22 




3 


PG 0244+194 


8.05 


-21.89 




3 


PG 0923+201 


8.95 


-22.37 


< 11.91 


3,7 


PG 0953+414 


8.40 


-21.98 




3 


PG 1012+008 


7.80 


-21.90 


< 11.48 


3,7 


PG 1029-140 


9.09 


-22.12 




4 


PG 1116+215 


8.22 


-22.82 


< 11.52 


4,5 


PG 1202+281 


8.30 


-21.92 


< 11.38 


4 


PG 1307+085 


7.86 


-22.02 


11.32 


4,6 


PG 1309+355 


8.01 


-22.12 


< 11.41 


4,5 


PG 1402+261 


7.30 


-20.92 


11.39 


4,5 


PG 1444+407 


8.08 


-22.02 


11.57 


4,5 


PG 1635+119 


8.11 


-22.17 




3 


Elliptical Galaxies 


NGC 821 


7.57 


-21.94 




8 


NGC 2778 


7.15 


-20.12 




8 


NGC 3377 


8.00 


-20.58 




8 


NGC 3608 


8.28 


-21.39 




8 


NGC 4291 


8.49 


-21.16 




8 


NGC 4473 


8.04 


-21.42 




8 


NGC 4564 


7.75 


-20.45 




8 


NGC 4649 


9.30 


-22.83 




8 


NGC 4697 


8.23 


-21.77 




8 


NGC 5845 


8.38 


-20.25 




8 


M 32 


6.59 


-17.36 




9 


NGC 3379 


8.00 


-21.47 




9 


NGC 4486B 


8.70 


-18.30 




9 


NGC 4742 


7.15 


-20.47 




9 


NGC 4261 


8.71 


-22.62 




9 


NGC 4374 


8.64 


-22.89 




9 


M87 


9.48 


-23.06 




9 


NGC 6251 


8.78 


-23.34 




9 


NGC 7052 


8.51 


-22.84 




9 


IC 1459 


8.30 


-22.92 




9 



Note. — Col.(l):source name; col. (2) BH mass; col. (3) i?-band absolute magnitude of host galaxies (mag); col. (4) far-infrared luminosity References — (1) 
Zheng ct al. 2002; (2) Vcilleux et al. 2002; (3) Dunlop ct al. 2003; (4) McLure fc Dunlop 2001; (5) Haas ct al. 2003; (6) Haas et al. 2000; (7) Sanders 
et al. 1989: (8) Gebhardt et al. 2003; (9) Kormcndy & Gebhardt 2003 
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Tabic 2. The Basic Physical Parameters of Type I ULIRGs 



IRAS 


z 


F A (5100(l + z)A) obs 


FWHM (Hfj) 


AL A (5100A) roBt 


F 2-io(,i+z)^/F 5imi+z) l 
(6) 


(1) 


(2) 


(3) 


(4) 


(5) 


IRAS F01572+0009 (Mrk 1014) 


0.163 


1.41 x 10" 1 " 


2100 


4.6 x 10 44 


5.3 x 10 3 


IRAS F07599+6508 


0.148 


1.50 x 10~ 15 


3150 


3.9 x 10 44 


9.3 x 10 1 


IRAS F11119+3257 


0.189 


1.45 x 10~ 16 


1980: 


6.7 x 10 43 


9.0 x 10 4 


IRAS Z11598-0112 


0.151 


1.72 x 10~ 16 


820 


4.7 x 10 43 


1.9 x 10 4 


IRAS F12540+5708 (Mrk 231) 


0.042 


1.25 x 10~ 14 


3130 


2.3 x 10 44 


3.4 x 10 2 


IRAS F13342+3932 


0.179 


1.44 x 10~ 16 


980 


5.6 x 10 43 




IRAS F15462-0450 


0.101 


7.40 x 10~ 17 


1460 


8.6 x 10 42 




IRAS F21219-1757 


0.113 


3.24 x 10~ 16 


2080 


8.0 x 10 43 





Note. — The prefix of the object name indicates the origin of IRAS fluxes. "F" refers to the IRAS Faint Source Catalogue, and "Z" means the Faint 
Source Reject File. The uncertainty on the FWHM is typically of order 10%; colons indicates values with a relative uncertainty of 20% (Zheng et al. 
2002). Col.(l):sourcc name; col.(2): rcdshift; col. (3):observed flux at 5100(1 + z)A [erg s" 1 cm- 2 A 1 ]; col. (4) FWHM of H/3 emission profile [km s ]; 
col. (5): continume luminosity [erg s~ 1 ] at rest frame; col. (6): Flux ratio F 2 — io(i+z)kcv/-^ , 5100 ( 1 _|_ ^2 — i0kcV(i+s) 1S unit of [erg s _1 cm -2 A 1 ] 
References. — (l)-(5): Zheng ct al. (2002) and their original spectra. (6) rGrorge et al. (2000) 



Table 3. The Parameters for the Extinction of Type I ULIRGs 



IRAS 


Af H [cm- 2 ] 


A v [mag] 


A« al " tlc [mag] 


references 


(1) 


(2) 


(3) 


(4) 


(5) 


IRAS F01572+0009 (Mrk 1014) 


< 2.6 x 10 2 " 


< 0.012 


0.087 


1 


IRAS F07599+6508 


> 10 24 




0.16 


2 


IRAS F11119+3257 


1.2 x 10 22 


0.57 


0.159 


1 


IRAS Z11598-0112 


6 x 10 20 


0.029 


0.066 


1 


IRAS F12540+5708 (Mrk 231) 


> 10 24 




0.078 


3 


IRAS F13342+3932 






0.012 




IRAS F15462-0450 






0.65 




IRAS F21219-1757 






0.21 





Note. — Col.(l):source name; col. (2): column density; col. (3): internal optical extinctions derived from column demsity. The value is evaluated by 
Ay/Nh — 4.8^3 4 6 1 X 10~ 23 magcm~ 2 (lcr dispersion) derived from AGN observation (Maiolino ct al. 2001b); col. (4) Galactic extinction at visual band 
References. — (1) Anabuki 2004; (2) Braito ct al. 2004; (3) Imanishi & Tcrashima 2004 



Table 4. The Basic Model Parameters (KUM03 Model) 



M gO [M ] 


fclGyr- 1 ] 


tw[yr] 


tcrit [yr] 


Across [yr] 


(1) 


(2) 


(3) 


(4) 


(5) 


10 1U 


14.7 


1.3 x 10 B 


4.5 x 10 8 


7 x 10 8 


10 11 


11.3 


3.5 x 10 s 


8 x 10 s 


1.2 x 10 9 


10 12 


8.6 


7 x 10 8 


1.2 x 10° 


1.8 x 10 9 


2 x 10 12 


8.1 


8 x 10 s 


1.4 x 10 9 


2 x 10° 



Note. — (1) M g o is the initial gas mass in galactic bulges. (2) k is the constant rate coefficient for SFR. (3) t w is the galactic wind time scale. (4) i cr it 
is the transition time from the host luminosity-dominat phase to the AGN luminosity dominant phase. (5) Across is the time when the BH growth stops. 
As for these parameters, see also text and the original paper (KUM03). 



